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BACKGROUND 

The  motions  of  a  cable  suspended  system  deployed  in  an  ocean  environment  has  long  been  of  interest 
to  the  oceanographic  community.  Programs  which  solve  for  the  steady  state  cable  configuration  for 
systems  deployed  in  a  variable  curient  profile  have  been  in  use  for  some  time  (reference  1,2).  Several 
programs  which  generate  the  relative  motions  of  suspended  bodies  as  a  result  of  a  random  surface  forcing 
function  are  also  in  existence.  In  general,  however,  these  programs  mod  a  I  only  one  specific  component  of 
the  deployed  system,  and  approximate  the  remainder  of  the  system.  A  program  which  could  generate  the 
dynamic  motion  of  an  arbitrary  cable  body  system  was  required. 


INTRODUCTION 


1  ?  CABUOY  is  a  Fortran  IV  programme veloped  for  the  Naval  Air  Development  Center  by  the  David  Taylor 
Naval  Ship  Research  and  Development  CenteTl^The  program  solves  for  the  two-dimensional  dynamic 
motions  of  a  cable  suspended  system  deployed  in  an  ocean  environmenvirefCT^cerSiVlntended  as  an 
engineering  tool,  the  program  can  model  any  combination  of  cable  segments  and  attached  bodies 
configured  as  either  a  towed,  moored,  or  free  drifting  system.  The  positions  of  all  bodies  and  the  resultant 
cable  tensions  are  calculated  as  functions  of  time.  This  time  domain  approach,  though  more  costly  in 
terms  of  execution  time  than  the  more  common  frequency  domain  method,  was  chosen  because  it  allows 
for  the  solution  of  transient  motions  which  may  induce  high  tension  loads  in  the  cable. 


V _ _^The  program  employs  a  finite  element  analysis  to  model  the  cable-body  system.  The  resultant 

differential  equations  are  integrated  using  a  Kutta-Merson  routine  which  automatically  reduces  the 
integration  time  step  until  specified  error  criteria  are  met-This  report  will  not  deal  further  with  the 
mathematics  involved  as  these  are  completely  described )i\  reference  (3).  Rather,  it  will  discuss  the 
program  input  parameters  and  how  they  are  manipulated  totfescribe  various  systems,  as  well  as  the 
information  available  in  the  output  data  listing.  \ 
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DISCUSSION 


SYSTEM  MODEL 

CA8U0Y  solves  for  the  dynamic  motions  of  moored,  drifting  or  towed  systems  consisting  of  a  surface 
float  and  any  combination  of  cable  segments  and  suspended  bodies  as  illustrated  in  figure  1 .  Free  drifting 
and  towed  cases  are  modeled  in  the  program  as  moored  systems  where  the  bottom  cable  segment  is 
defined  to  be  massless  and  highly  elastic.  The  towed  system  is  further  considered  as  a  free  floating  system 
where  the  top  of  the  cable  is  restricted  from  moving  horizontally.  The  system  is  then  "towed"  by  inputina  a 
coplanar  current  velocity  of  magnitude  equal  to  the  tow  velocity.  The  tensions  measured  in  the  fictitious 
cable  segment  of  the  drifting  system  modal  are  generally  less  than  1  percent  of  the  total  system  weight 
and  have  negligible  effects  on  the  dynamic  solution. 

CABUOY  allows  the  user  to  arbitrarily  define  the  system's  initial  position  and  velocity  from  which  the 
dynamic  data  is  calculated.  Thus,  it  is  possible  to  compute  the  time  required  for  a  system  to  reach  an 
equilibrium  condition  from  some  initial  configuration.  Alternatively,  the  initial  conditions  can  be 
automatically  set  to  the  steady-state  values  calculated  by  CABUOY.  This  minimizes  system  transients  and 
is  useful  for  generating  long  term  dynamic  data.  Such  versatility  greatly  increases  the  usefulness  of  the 
program. 

A.  Surface  Float 

Two  classes  of  surface  floats  can  be  modeled  by  CABUOY:  spar  buoys,  defined  as  buoys  with  large 
draft  to  diameter  ratios,  and  spheroidal  buoys  whose  limiting  cases  range  from  a  thin  disc  to  a  spar.  These 
two  classes  cover  the  spectrum  of  most  commonly  used  surface  floats. 

During  the  dynamic  calculation  phase,  CABUOY  generates  the  coefficients  for  pitch,  heave,  and  surge, 
and  then  calculates  the  response  of  the  float  to  the  surface  wave  forcing  function.  Due  to  the  large 
amounts  of  computer  time  normally  involved  in  these  calculations,  several  alternative  approximations 
have  been  included  which  reduce  execution  time  at  the  expense  of  accuracy.  When  using  these 
approximations,  the  surface  float  drag  areas  and  added  mass  remain  constant  at  the  steadv-state  vaiues  or 
are  calculated  only  at  the  end  of  each  print  interval.  These  approximations,  though  not  accurate,  are  useful 
when  evaluating  the  effects  of  subsurface  assembly  changes  on  system  performance  where  many  runs 
may  be  necessary. 

The  user  also  has  the  option  of  coupling  the  surface  float  directly  to  the  ocean  surface,  eliminating  the 
interaction  between  waves  and  float.  This  mode  is  useful  in  representing  a  system  deployed  off  of  a  ship 
when  the  ship  motion  is  known.  The  typical  laboratory  arrangement  where  an  oscillating  arm  drives  the 
top  of  the  cable  and  the  response  of  the  lower  unit  is  measured  can  also  be  simulated  in  this  manner. 


B.  Cable 

In  CABUOY,  the  system  suspension  is  divided  into  a  number  of  rigid  extensible  segments  whose  end 
points  are  defined  by  the  bodies  suspended  on  the  cab'e.  To  minimize  the  errors  induced  by  the 
assumption  of  a  rigid  cable,  long  continuous  cable  length  can  be  divided  up  into  several  segments.  This  is 
especially  useful  at  the  top  of  the  system  where  the  cable  catenary  is  greatest. 
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MODEL  DEFINITION 


Figure  1  -  CABUOY  {Model  Definition) 
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The  stress-strain  relationship  for  each  cable  segment  is  defined  by: 

T  =  C1*f  C2  +  CINT*c  +  TREF 
where  T  =  cable  tension 

Cj  C2_  are  constants 
t  =  strain 

CINT  =  Internal  Damping  Coefficient 
*  =  time  rate  of  change  of  strain 
TREF  =  constant 


This  equation  allows  CABUOY  to  model  nonlinear  cables,  such  as  bungee  cords,  as  well  as  linear 
materials  such  as  steel  and  Kevlar,  through  the  appropriate  choice  of  coefficients. 

C.  Suspended  Bodies 

All  suspended  bodies  are  considered  to  be  point  masses  of  known  drag  area  and  virtual  mass,  with  all 
forces  acting  at  the  center  of  gravity  of  the  body  Because  of  their  wide  use  as  vertical  motion  isolation  in 
sonobuoy  systems,  a  special  routine  has  been  included  for  flat  discs  which  calculates  the  virtual  mass  as  a 
function  of  the  oscillation  amplitude.  For  cases  where  several  bodies  are  connected  by  short  cable 
segments,  such  as  in  an  instrumented  line  array,  it  is  advantageous  to  input  the  bodies  as  a  single  element 
with  a  drag  area  equal  to  that  of  the  entire  array.  This  minimizes  the  number  of  cable  segments  and  thus 
the  execution  time  without  significantly  increasing  the  computational  error. 


ENVIRONMENT 

The  ocean  environment  utilized  in  CABUOY  consists  of  a  user-defined  sea  surface  which  decays  ex¬ 
ponentially  with  depth  in  conjunction  with  a  variable  current  velocitv  profile.  The  ocean  surface  is 
generated  by  the  summation  of  N  sine  waves,  the  amplitudes,  frequencies,  and  phases  of  which  may  be 
defined  by  the  user  or  by  an  internal  wave  generator.  The  current  profile  is  input  as  a  table  of  velocity 
versus  depth.  For  the  moored  system,  the  depth  of  the  last  current  velocity  is  defined  to  be  the  ocean 
bottom.  For  a  free  floating  system,  the  current  velocity  at  all  depths  greater  than  the  maximum  depth 
specified  is  equal  to  the  last  current  velocity  input. 


DATA  INPUT 

CABUOY  input  data  is  arranged  as  illustrated  in  figure  2.  Program  variables  are  defined  in  appendix  A, 
along  with  a  data  format  guide.  In  general,  data  input  is  a  straight  forward  process.  However,  because  of 
the  various  systems  CABUOY  is  designed  to  evaluate,  several  variables  have  multiple  definitions, 
depending  on  which  system  is  being  analyzed.  This  section  will  review  just  those  variables. 
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Figure  2  -  CABUOV  (Data  Input) 
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A.  FSM 

The  motion  at  the  top  of  the  cable  is  dependent  on  the  variable  FSM  (figure  3).  When  FSM  is  less  than 
1C00,  the  motion  is  defined  by  a  summation  of  sine  waves  with  a  frequency,  amplitude,  and  phase 
determined  by  the  variables  AXSM,  AYSM,  FSM,  and  FIDSM.  This  motion  acts  directly  at  the  top  of  the 
cable  system  and  eliminates  the  surface  float  from  the  calculations.  For  1 000  <  FSM  <  2000.  the  motion 
at  the  top  of  the  cable  is  set  equal  to  the  surface  waves  defined  by  FSW,  ASW,  and  FIDSW.  The  surface 
float  is  again  eliminated  from  the  calculations. 

In  those  cases  where  2000<  FSM  <  2900,  the  program  considers  the  surface  float  to  be  a  spar  buoy, 
and  calculates  the  response  of  the  cable  buoy  system  to  the  surface  wave  forcing  function.  The  float  is 
defined  by  considering  AXSM  (k)  to  be  the  cross  sectional  area  of  the  buoy  at  a  vertical  distance  AYSM  (k) 
from  the  origin  of  the  local  buoy  coordinate  system.  All  inertia  and  coupling  coefficients  are  calculated  by 
the  program.  The  origin  of  the  float  coordinate  system  can  be  anywhere  the  user  defines,  and  AYSM 
should  be  negative  or  positive  as  required  to  define  areas  above  and  below  the  origin. 

Where  2900<i  FSM  <  3000,  the  buoy  is  considered  to  be  an  arbitrary  shape  defined  in  the  same 
manner  as  a  spar  buoy.  However,  all  inertia  and  coupling  coefficients  are  defined  by  the  user.  For  3000  < 
FMS  <  4000,  the  program  considers  the  float  tc  oe  a  spheroid.  AXSM  (1 )  and  AYSM  (1 )  are  considered  to 
be  the  horizontal  and  vertical  semi-axes  respectfully.  All  inertia  and  coupling  coefficients  are  again 
calculated  by  the  program.  The  origin  is  considered  to  be  at  the  geometric  center  of  the  float. 

Additional  data  required  to  define  the  surface  float  as  well  as  any  subsurface  package  are  input  on 
separate  cards.  All  required  variables  are  defined  in  appendix  A.  Figure  4  illustrates  the  various  forces  and 
their  geometric  relationship  to  the  float  Cg  It  should  be  noted  that  the  ir  put  variable  TIY  represents  the 
force  at  the  top  of  the  cat  a.  As  this  force  is  in  a  vertical  up  direction,  the  va  ue  enter  for  TIY  is  the  negative 
of  the  actual  value. 

B.  ASW 

The  surface  wave  forcing  function  is  defined  by  ASW  (figure  5).  When  ASW  <  1000,  the  user  inputs 
the  amplitude,  frequency,  and  phase  of  up  to  20  smusodial  wave  components  which  are  then  summed 
together  to  generate  the  surface  wave.  For  1 000  <  ASW  <  2000,  the  program  generates  a  random  sea 
surface  with  wave  components  calculated  fiom  the  Pierson-MosKowitz  energy  spectrum  for  a  well 
developed  sea.  ASW  (1 )  - 1000  is  equal  to  the  significant  wave  height,  and  FRSW  (1  )and  FRSW  (2)  are  the 
lower  and  upper  frequency  limits  for  the  spectrum.  All  components  are  phase  shifted  relative  to  each  other 
by  360/NSW  degrees. 

For  the  case  where  2003  <  ASW  <  2008,  a  standard  random  sea  surface  is  generated  using  wave 
components  defined  by  the  program  reference  (4).  ASW  -  2000  is  equal  to  the  desired  sea  state. 


C.  Output  Time  Increments 

The  variables  DTI .  DT2.  TINV1 ,  and  T0TT  are  used  to  control  the  output  of  the  dynamic  solution 
(figure  6).  TINV1  defines  an  initial  interval  during  which  data  is  printed  every  DTI  seconds.  This  interval 
could  be  used  to  give  limited  information  during  a  period  where  transient  system  motions  are  decaying 
For  all  times  after  T.INV1  until  TOTT,  data  is  printed  out  every  DT2  seconds.  This  would  be  wnere  fine 
resolution  is  required  to  analyze  long  term  motions.  The  time  increments  could  be  changed  if  the  user  was 
more  interested  in  the  initial  system  response,  and  how  much  time  was  required  to  reach  steady  state 
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Figure  3  -  CABUOY  (Surface  Motion) 
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Figure  6  -  C/._JOY  (Time  Control) 
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D.  Cable  Stress  Versus  Strain 

CABUOY  uses  the  generalized  equation: 

Tension  =  TREF  +  C1*cc2  +  CINT  *  * 

to  relate  cable  stress  and  strain  (figure  7).  This  equation  allows  the  user  to  closely  model  the  behavior  of 
many  varied  cable  materials.  Fora  linearly  elastic  material,  C2  =  1  and  Cl  =  AE..where  A  is  the  cable  cross 
sectional  area  and  E  is  the  elastic  modulus.  Fora  nonlinear  material,  the  user  can  select  values  of  Cl  and 
C2  to  closely  model  the  material  in  the  immediate  area  of  interest.  It  should  be  noted  that  if  TREF  is  non¬ 
zero,  the  cable  length  input  should  be  at  a  tension  of  TREF.  TMIN  must  also  be  defined  to  prevent  the  cable 
from  being  compressed  during  the  dynamic  phase. 


E.  Bodies 

CABUOY  considers  each  body  to  be  a  point  mass  with  a  fixed  drag  area.  For  those  special  cases  where 
the  body  is  a  circular  damping  disc,  a  special  routine  is  available  which  calculates  the  effective  drag  and 
added  mass  terms  as  a  function  of  the  oscillation  amplitude.  The  user  enables  this  routine  by  entering  the 
negative  of  the  physical  drag  area  for  the  variable  CDABX  if  the  disc  is  oriented  in  the  X-Z  plane,  or  the 
variable  CDABY  if  the  disc  is  in  the  Y-Z  plane.  The  remaining  drag  area  would  be  entered  normally. 

DATA  OUTPUT 

Data  output  is  formatted  as  illustrated  in  figure  8.  First,  all  input  data  are  listed.  The  program  calculates  and 
outputs  the  surface  float  constants  required  for  the  dynamic  calculations  based  upon  the  draft  estimate 
from  the  input  data.  Using  these  values  as  initial  approximations,  the  program  calculates  the  steady-state 
configuration  of  the  system.  These  steady-state  calculations  use  the  finite  element  analysis  techniques 
and  iteration  routines  from  the  programs  FF2E  and  MR3E.  and  considers  the  cables  to  be  flexible 
extensible  segments.  The  results  are  more  accurate  than  would  be  obtainable  using  the  rigid  cable 
segments.  The  results  from  each  iteration  are  printed  out,  and  when  a  solution  is  reached  the  constants  for 
the  surface  buoy  are  recomputed  and  output  for  the  actual  float  draft.  The  exact  surface  float  tilt  is  written, 
as  well  as  the  approximate  float  tilt  used  in  the  dynamic  calculations.  The  approximate  tilt  may  vary  from 
the  exact  tilt  as  a  result  of  the  cable  assumptions  (rigid  extensible  segment)  made  for  the  dynamic  phase. 

The  program  then  lists  the  steady-state  configuration  of  the  cable  system.  This  represents  the  "exact" 
steady-state  system  configuration.  By  fitting  the  rigid  CABUOY  cable  segments  to  the  exact  solution,  the 
program  generates  a  set  of  initial  cable  conditions  which  will  be  used  in  the  dynamic  calculations  if 
required  by  the  user.  The  program  then  calculates  the  cable  motions  and  tensions  and  outputs  these  data 
as  a  function  of  time.  Data  are  output  for  the  surface  wave  (WAVE),  surface  float  center  of  gravity  (BUOY), 
top  of  the  cable  (Node  0),  and  the  bottom  of  each  cable  segment  (Nodes  1 ,2,3  etc).  No  data  are  output  for 
the  bottom  of  the  last  cable  segment,  which  is  assumed  to  be  fixed  in  space.  While  the  program  computes 
the  dynamic  motions,  it  also  stores  the  percentage  of  the  time  that  float  draft  and  pitch  exceed  predefined 
limits.  These  percentages  are  then  output  at  each  print  interval. 

CABUOY  EXAMPLE 

To  demonstrate  the  application  of  CABUOY  for  the  solution  of  a  real  world  problem,  the  moored 
system  illustrated  in  figure  9  was  evaluated  and  the  input/output  data  listed  in  tables  I  and  II.  It  is 
suggested  that  the  example  be  closely  studied  and  understood  before  the  user  attempts  a  problem  of  his 
own. 
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Figure  8  -  CABUOY  (Data  Output) 
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Speroidal  Surface  Buoy 
Waterplane  radius 
Ora  ft 

Weight  m  air 
Moment  of  inertia 
Submerged  x-drag  area 
Submerged  y*drag  area 


Wind  Load  4  5  lb  (20  03N) 


YCG  =  3.7  ft  (1  13  m) 


Fluid  density  =  1  99  slugs/ft^  (1026  kg/m^) 


0.53  ft  (0.162  m) 

5  3  ft  (1  62m) 

38  lb  069.1  N) 

1  8  slug  ft^  (2  44  kg  m^) 
3  15  ft  2  (0  293  m2) 

0  22  ft2  (0.020  m2) 


Body  1  Body  2  Body  3 


Weight  in  fluid,  lb  (4  45N) 

-40 

0 

100 

1 5.C 

x-drag  area,  ft2  (0  0929  m2) 

0. 

44 

0  18 

0.70 

y-drag  area,  ft2  (0.0929  m2) 

0 

44 

1.95 

0  08 

x-virtual  mass,  slugs  (14  6  kg) 

0 

73 

0.55 

1  92 

y-vutual  mass,  slugs  (14  6  kg) 

0 

73 

2.75 

1  35 

Cable  2 

: 7777777777777 


Cable  1 

Cable  2 

Diameter,  in.  (2.54  cm) 

0  15 

0  25 

Normal  drag  coefficient 

1  4 

1  4 

Tangential  drag  coefficient 

0  02 

0  02 

Weight  m  fluid.  Ib/ft  (14.6  N/m) 

004 

0  02 

Mass,  slugs/ft  (47  84  kg/m) 

00015 

0  0013 

Reference  tension.  Ib  (4  45  N) 

0 

100. 

C,  Ib  (4  45  N) 

22.000. 

1.600 

c2 

1  0 

1  5 

Internal  damping  coefficient. 

0. 

0 

Ib  sec  (4  45  N  sec) 

Profile-  Depth 

Current 

(ft) 

(knots) 

(0.305  m) 

(0  51 5  m/s) 

0 

2.5C 

500 

1  30 

1000. 

0  50 

1500. 

0  50 

Figure  9  -  Parameters  for  Moored  Buoy-Cable- Body  System 
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Program  Time  Requirement 

As  noted  in  the  original  NSRDC  report,  CABUOY  requires  a  large  amount  of  memory  and  may  take 
several  minutes  of  CPU  time  to  generate  a  solution.  Generally,  the  number  of  surface  waves  or  motions 
(NSW  or  NSM)  and  cable  segments  (NCAB),  the  stiffness  of  the  cable  segments  (AE),  the  degree  of 
surface  float  approximations  (IBUOY),  and  the  amount  of  output  data  required  (TOTT)  are  the  primary 
factors  affecting  the  execution  time.  The  following  techniques  will  reduce  the  required  computer  time. 

NSW/NSM  -  Reducing  the  number  of  surface  wave  components  will  reduce  execution  time.  Generally, 
it  is  possible  to  evaluate  systems  response  to  a  simple  harmonic  motion.  Only  in  selected  cases  would  the 
user  evaluate  system  response  to  a  random  sea  surface. 

NCAB  -  Reducing  the  number  of  cable  segments  will  directly  reduce  execution  times  (NCAB  >  2). 
Generally,  nearsurface  cable  segments  should  be  divided  into  small  lengths,  while  deep  segments  can  be 
grouped  into  one  long  segment.  Often  a  vertical  hydrophone  array  can  be  approximated  by  one  cable 
segment  where  cable  drag  area  is  defined  to  be  equivalent  to  the  array  drag,  with  appreciable  reductions 
in  CPU  time. 

AE  -  Stiff  cable  segments  (high  equivalent  spring  constants)  require  longer  execution  time  because  of 
the  smaller  time  increments  used  during  the  integrations.  If  possible,  it  would  be  more  efficient  to  input  a 
shorter,  more  elastic  segment  which  will  stretch  to  the  original  cable  segment  length.  The  effects  of  this 
change  on  the  system  response  will  have  to  be  evaluated  by  the  user. 

IBUOY  -  Generally,  the  more  often  surface  float  constants  have  to  be  calculated,  the  more  execution 
time  is  required.  Therefore,  unless  the  accuracy  is  required,  the  user  can  economize  by  setting  IBUOY 
equal  to  0  or  1 .  Setting  IBUOY  equal  to  -1 0  should  be  avoided  except  in  cases  where  surface  float  motion 
is  limited. 


TOTT  -  Obviously,  the  more  output  data  required,  the  longer  the  program  will  run.  Execution  is  not 
affected  by  increasing  DTI  or  DT2.  as  the  integration  time  step  is  determined  by  the  computer,  and  this 
determines  execution  time.  The  only  charges  which  change  with  DTI  and  DT2  are  those  associated  with 
line  printers. 
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TABLE  I 

INPUT  PARAMETERS 

Col  FMT  Variable  Data 

1-3  13  NCASES  001 

1-80  20A4  TITLE  CABUOY-CASE  2  of  NSRDC  Report 

1-3  13  NSM  001 

4-8  13  NSW  001 

7-9  13  NCAB  005 

10- 12  13  NCUR  004 

13-15  13  ITER  002 

16-18  13  MTRC  000 

19-21  13  I  BUOY  -01 

22-24  13  ISPAR  0 

1-10  FI  0.4  AXSM(I)  .530 

11- 20  FI  0.4  AYSM  (I)  5.3 

21-30  FI  0.4  FSM  (I)  3500. 

31-40  FI 0.4  FIDSM  (I)  1.0 

I- 10  FI  0.4  ASW(J)  7.5 

II- 20  FI  0.4  FRSW  (J)  .10 

21-30  FI  0.4  FIDSW(J)  0.0 

I- 10  FI  0.4  RHO  1.99 

II- 20  FI  0.4  AMC  1.0 

21-30  FI  0.4  AFAC  1.0 

31-40  FI  0.4  TBH  .0 

41-50  FI  0.4  TBYMX  999999. 

I- 10  FI  0.4  SUBM  0.0 

II- 20  FI  0.4  TWX  4.5 

21-30  FI  0.4  TIY  -100. 

31-40  FI  0.4  TMIN  0. 

I- 10  FI  0.4  CDASPX  0.0 

II- 20  FI  0.4  CDASBX  6.3 

21-30  FI  0.4  VSP  0.0 

31-40  FI  0.4  SPXK  0.0 

41-50  FI  0.4  SPYK  0.0 

51-60  FI  0.4  YYSP  0.0 

I- 10  F10.4  TINVI  10.0 

II- 20  FI  0.4  DTI  0.1 

21-30  FI  0.4  TOTT  50.0 

31-40  FI  0.4  DT2  0.5 


1-0 

FI  0.2 

FLC  (K) 

200. 

700 

900. 

500. 

200. 

11-20 

FI  0.4 

DCI  (K) 

.15 

.15 

.15 

.25 

.25 

21-30 

FI  0.4 

CDN  (K) 

1.4 

1.4 

1.4 

1.4 

1.4 

31-40 

FI  0.4 

CDT  (K) 

.02 

.02 

.02 

.02 

.02 

41-50 

FI  0.4 

WC  (K) 

.04 

.04 

.04 

.02 

.02 

51-60 

FI  0.6 

CM  (K) 

.0015 

.0015 

.0015 

.0013 

.0013 

1-10 

FI  0.0 

Cl  (K) 

22000. 

22000. 

22000. 

1600. 

1600. 

11-20 

FI  0.4 

C2  (K) 

1.0 

1.0 

1.0 

1.5 

1.5 

21-30 

F10.4 

CINT  (K) 

0. 

0. 

0. 

0. 

0. 

31-40 

FI  0.2 

TREF  (K) 

0. 

0. 

0. 

100. 

100. 

41-50 

F10.2 

TENI (K) 

0. 

0. 

0. 

0. 

0. 

51-60 

FI  0.4 

PHID  (K) 

9999. 

0. 

0. 

0. 

0. 

61-70 

FI  0.4 

XPI  (K) 

0. 

0. 

0. 

0. 

0. 

71-80 

FI  0.4 

YPI  (K) 

0. 

0. 

0. 

0. 

0. 
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APPENDIX  A 

DEFINITION  OF  INPUT  VARIABLES 


NCASES 

TITLE 

NSM 


NSW 

NCAB 

NCUR 

ITER 


MTRC 

IBUOY 


ISPAR 


AXSM 

AYSM 

FSM 

FIDSM 


—  Number  of  cases  to  be  run.  One  case  is  defined  by  data  cards  2  through  17. 
NCASES  > 1 

—  Title  of  data  case.  Title  is  printed  on  output  for  identification. 

—  Number  of  surface  motion  components.  For  1000  <  FMS  (1)  <  2000 

Surface  motion  =  Surface  Wave,  NSM  is  automatically  set  equal  to  NSW.  1  <  NSM 
<  20 

—  Number  of  surface  wave  components.  1  <  NSW<  20 

—  Number  of  cable  segments.  2  <  NCAB  ^  50 

—  Number  of  current  profile  data  points.  2  £  NCUR  <1 10 

—  ITER  =  0  prescribed  initial  steady  state  conditions 

=  1  free  floating  cable  system 
=  2  moored  cable  with  given  length  in  given  depth 
=  3  towed  cable  system 

=  4  iteration  scheme  to  be  programmed  by  user 

—  MTRC  <  0  if  input  data  are  entered  in  English  units. 

>  1  if  input  data  are  entered  in  metric  units. 

—  IBUOY  <-10  Surface  buoy  drag  areas  and  added  masses  remain  at  constant  values,  and 

displaced  volume  is  equal  to  the  draft  times  the  steady-state  waterplane 
area. 

<  -  1  Surface  buoy  drag  areas  and  added  masses  remain  at  constant  values 

corresponding  to  the  steady-stae  draft. 

=  0  Surface  buoy  drag  areas  and  added  masses  are  updated  at  the  end  of  each 
print  interval. 

<  1  Surface  buoy  drag  areas  and  added  masses  are  updated  continuously  as  a 

function  of  buoy  submergence. 

Degree  of  polynomial  used  to  approximate  variation  of  cross  sectional  area  of  a  spar 
buov.  ISPAR  <  NSM-1 

FSM  (1  j  <  1000  Amplitude,  frequency  and  phase  of  surface  motion  which  is 
defined  by  the  equations: 

NSM 

xsm  =  I  AXSM  (k)  *  COS  (-2X*FSM  (k)  *t  +  FIDSM  (k)  *x/180) 
k=1 


NSM 

Ysm  =  I  -AYSM  (K)  *SIN  (-2x*FSM  (k)  *  t  +  FIDSM  (k)  *x/180) 
k=1 

(ft.,  hz  8-  deg  or  m,  hz  &  deg) 

1 000  <  FSM  (1 )  <  2000  =  Prescribed  surface  motion  components  are  set  equal  to 
the  surface  wave  components  by  setting  AXSM  (k)  =  AYSM  (k)  =  ASW  (k), 
FSM  (k)  =  FRSW  (k),  FIDSM  (k)  =  FIDSW  (k)  and  NSM  =  NSW. 


2000  <  FSM  (1 )  <  2900  -  Surface  buoy  is  considered  to  be  a  spar  buoy.  Program 
calculates  the  inertia  and  forcing  coefficients.  AXSM  (k)  is  the  cross- 
sectional  area  of  the  buoy  at  a  vertical  distance  AYSM  (k)  from  the  origin  of 
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ASW 

FRSW 

FIDSW 


RHO 

AMC 

AFAC 

TBH 

TBYMX 

SUBM 

TWX 

TIY 


APPENDIX  A  (Cont'd) 


the  local  buoy  coordinate  system.  AYSM  (NSM)-AYSM  (1 )  is  the  total  length 
of  the  buoy.  FIDSM(1)is  the  added  mass  coefficient  for  surge,  Ks.  FIDSM(2) 
is  the  ratio  a  =  Isw/  (xjr^/4),  where  Isw  is  the  moment  of  inertia  of  the  water 
plan  area  and  r^ is  the  radius  of  the  water  plane  cross  section.  For  a  circular 
section  spar  buoy,  Ks  =  a  =  1 .  These  values  are  automatically  preset  if 
FI  DSM  (1 )  =  FIDr<  M  (1 )  =  0.  All  inertia  and  forcing  coefficients  are  calculated 
by  the  program. 

2900  <  FSM  (1 )  <  3000  Surface  buoy  is  considered  to  be  an  arbitrary  shape.  Inertia 
and  forcing  coefficients  are  input  as  data,  AXSM  (k),  AYSM  (k),  FSM  (k)  and 
FIDSM  (k)  are  defined  as  for  a  spar  buoy. 

3000  <  FSM  <  4000  Surface  buoy  is  considered  to  be  a  spheroid.  AXSM  (1 )  and 
AYSM  (1)  are  the  horizontal  and  vertical  semi-axes.  All  other  variables 
should  be  set  equal  to  0. 

ASW  <  1 000.  amplitude,  frequency  and  phase  of  surface  wave  components,  which 
are  defined  by  the  equation: 

v  NSW 

xsw  =  r  ASW  (k)  *  COS  (-2x*FRSW  (k)  *t  +  FIDSW  (k)  *x/180). 
k=1 
NSW 

Ysw  =  r  -ASW  (k)  *SIN  (-2X* FRSW  (k)  *t  +  FIDSW  (k)  *A/180). 
k=1 

(ft.  Hz  &  deg  or  m,  Hz  &  deg). 

1 000  ^  ASW  (1 )  <  2000  sea  surface  is  generated  using  the  Pierson-Moskowitz  sea 
spectrum.  ASW  (1 )  -1 000  is  the  significant  wave  height,  and  FRSW  (1 )  and 
FRSW  (2)  are  the  lower  and  upper  frequencies  of  the  spectrum.  Wave 
components  are  phase  shifted  by  360/NSW  degrees,  with  FIDSW  (1 )  being 
the  phase  of  the  lowest  frequency  component. 

2003  <  ASW  (1 )  <  2008  A  standardized  random  sea  surface  is  generated  using 
internally  generated  wave  components. 

ASW  (1 )  -2000.  is  the  desired  sea  state.  Sea  states  3  through  8  are  available. 
Fluid  density  (Slugs/ft3  or  kg/m3j 

Added  mass  coefficient  for  cable.  AMC  =1.0  for  round  cable 
Cross-sectional  correction  factor.  AFAC  =  1 .0  for  cylindrical  cable.  AFAC  =  cable 
cross-sectional  area  /  (*d^/4)  for  non-round  cable. 

Applied  force  on  lower  weight  in  x-  direction  (lbs  or  N) 

Maximum  absolute  value  of  tension  in  cable  just  below  buoy.  For  buoy-cable 
systems,  set  TBYMX  =  99999.  (lbs  or  N) 

Depth  at  which  current  velocity  is  measured  for  surface  float  drag  calculations 
(ft.  or  m). 

Horizontal  force  acting  at  the  top  of  the  cable  (lb  or  N). 

For  cases  with  a  float,  TWX  represents  the  wind  loading. 

Vertical  component  of  tension  at  top  of  cable  (lb  or  N). 
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TMIN 

CDASPX 

CDAS8X 

VSP 

SPXK 

SPYK 

YYSP 


TINV1 

DTI 

TOTT 

DT2 

FLC 

DCI 

CDN 

CDT 

WC 

CM 

Cl 

C2 

CINT 


XPI 

YPI 

WBD 

CDABX 

CDABY 

XMBV 

YMBV 


For  cases  with  a  float,  T!Y  =-(steady-state  buoyancy  of  float). 

Minimum  tension  which  can  be  supported  by  cable,  (lb  or  N). 

Drag  area  of  surface  package  in  x  direction  (ft^  or  M^). 

Drag  area  of  surface  buoy,  excluding  surface  package,  in  x  direction  (ft^  or  M^). 
Volume  of  surface  package  (ft  ^  or  M^). 

Surface  package  added  mass  coefficient  for  surge. 

Surface  package  added  mass  coefficient  for  heave. 

Vertical  distance  of  center  of  buoyancy  of  surface  package  measured  from  the  origin 
of  the  local  buoy  coordinates  (ft  or  M). 

Initial  print  time  interval  (sec). 

Print  time  step  during  initial  interval  (sec). 

Total  time  for  which  dynamic  calculations  are  desired  (sec). 

Print  time  step  during  final  interval  (sec). 

Length  of  cable  segment  (ft  or  m). 

Diameter  of  cable  segment  (in.  or  cm). 

Normal  drag  coefficient  of  cable  segment. 

Tangential  drag  coefficient  of  cable  segment. 

Weight  in  fluid  of  the  cable  segment  at  the  reference  cable  tension  (*^s/ft  or  n/m). 
Mass  of  cable  segment  at  the  reference  cable  tension  (slugs/ft  or  kg/m). 

The  stress-strain  curve  for  the  segment  is  defined  by  the 
equation 

Tension  =  TREF  4-  Cl  *cc2+  CINT  *  * 

For  linearly  elastic  materials.  Cl  =  AE  and  C2  =  1. 

For  free  floating  and  towed  cases  (ITER  =  1  or  3)  the  last  cable  segment  connecting 
the  lower  weight  to  the  ocean  bottom  is  fictitious.  In  these  cases.  Cl  =  DCI  =  CDN  = 
CDT  =  WC  =  CM  =  CINT  =  0,  FLC  =  2  *FLC  (NCAB  -  1 )  and  C2  =  1. 

These  values  are  set  automatically.  Blank  cards  should  be  inserted  for  data. 

PHID  <  360  TENI  is  the  initial  value  of  cable  tension  (lbs  or  N). 

PHID  is  initial  value  of  cable  tilt  (<J>)  (deg). 

PHID  >  360  Initial  values  of  tension  and  cable  tilt  are  set  equal  to  steady-state 
values  calculated  by  the  program.  In  this  case,  0  can  be  entered  for  TENI. 

Initial  value  of  x  (ft/sec  or  m/sec). 

Initial  value  of  y  (ft/sec  or  m/sec). 

Weight  in  fluid  of  body  (lb  or  N). 

Drag  area  of  body  in  x  direction  (ft2  or  m2). 

Drag  area  of  body  in  y  direction  (ft^  or  m2). 

Virtual  mass  of  body  in  x  direction  (slugs  or  kg). 

Virtual  mass  of  body  in  Y  direction  (slugs  or  kg). 

If  CDABX  is  negative,  the  program  considers  the  body  to  be  a  circular  disk  with  plane 
perpendicular  to  the  x-axis,  and  calculates  the  drag  and  added  mass,  using 
I  CDABX  I  as  the  actual  drag  area  and  XMBV  as  the  actual  mass  of  the  disk. 
CDABY  and  YMBV  retain  their  original  definitions.  Similar  remarks  apply  if  CDABY  is 
negative,  except  that  the  plane  of  the  disk  is  now  perpendicular  to  the  Y  axis. 
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YY 

CCK 

CDASY 

WAST 

RWY 

RTX 

RTY 

YCG 

BINT 
XSI  I 
ZETI  > 
SYDl) 


XPSI  ) 
ZTPI  } 
SYPDI  ) 
DFTLIM 
SYDLIM 


Depth  of  current  velocity  data  point  (ft  or  m). 

For  moored  case,  YY  (NCUR)  is  the  ocean  bottom. 

Velocity  of  current  at  depth  YY  (kts  or  m/sec). 

Drag  drea  of  surface  buoy  in  Y  direction  (ft^  or  m^). 

Weight  in  air  of  buoy  and  surface  package  (lb  or  N). 

Vertical  distance  of  wind  loading  center  of  pressure  from  buoy  center  of  gravity 
(YCG)  (ft  or  m). 

Distance  of  cable  attachment  point  from  YCG  in  x  direction  (ft  or  m). 

Distance  of  cable  attachment  point  from  YCG  in  y  direction  (ft  or  m). 

Vertical  distance  of  center  of  gravity  of  buoy  and  surface  package  measured  from  the 
origin  of  the  local  buoy  coordinate  system  (ft  or  m). 

Moment  of  inertia  of  buoy  and  surface  package  about  YCG  (slug  ft  ^  or  km^). 
Initial  values  of  X.<,4/)>  where  <  is  the  vertical  displacement  of  the  float  center 
of  gravity  from  its  equilibrium  position  (ft,  ft  &  deg  or  m,  m  &  deg). 

For  SYDI_>  360,  the  program  sets  the  initial  values  for  buoy  inclination.draft  and 
vertical  velocity  equal  to  the  steady-state  values  previously  calculated  by  the 
program.  This  will  minimize  transient  dynamic  motions  of  the  surface  body. 
Initial  values  of  (iXM-  (ft/sec,  ft/sec.  Deg/sec  or  m/sec,  m/sec,  deg/sec). 


Limiting  value  of  draft  and  pitch  for  which  the  program  calculates  the  percent  of 
time  these  values  are  exceeded. 


(The  following  variables  are  described  in  NSRDC/SPD-0633-02) 


AKZ  Added  mass  coefficient  in  heave  for  arbitrary  buoy. 

AXP  A^/p  VL  for  arbitrary  buoy 

APP  (A^-BINTJ/pVL^  for  arbitrary  buoy 

AFKX  FKx/pV-  for  arbitary  buoy 

i*tu 

*FKZ  FKyA»Vyw  for  arbitrary  buoy 

AFKP 
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CABUOY 

INPUT  PARAMETERS 
DATA  FORMAT  GUIDE 


Card 

1 

2 

3 


4 


5 

6 


7 


8 


9 


10 


Col 

FMT 

1-3 

13 

1-80 

20A4 

1-3 

13 

4-6 

13 

7-9 

13 

10-12 

13 

13-15 

13 

16-18 

13 

19-21 

13 

22-24 

13 

1-10 

FI  0.4 

11-20 

FI  0.4 

21-30 

FI  0.4 

31-40 

FI  0.4 

1-10 

FI  0.4 

11-20 

FI  0.4 

21-30 

FI  0.4 

1-10 

FI  0.4 

11-20 

FI  0.4 

21-30 

FI  0.4 

31-40 

FI  0.4 

41-50 

FI  0.4 

1-10 

FI  0.4 

11-20 

FI  0.4 

21-30 

FI  0.4 

31-40 

FI  0.4 

1-10 

FI  0.4 

11-20 

FI  0.4 

21-30 

FI  0.4 

31-40 

FI  0.4 

41-50 

FI  0  4 

51-60 

FI  0.4 

1-10 

FI  0.4 

11-20 

FI  0.4 

21-30 

FI  0.4 

31-40 

FI  0.4 

1-10 

FI  0.2 

11-20 

FI  0.4 

21-30 

FI  0.4 

31-40 

FI  0.4 

41-50 

FI  0.4 

51-60 

FI  0.6 

B-1 


Variable 

N CASES 

TITLE 

NSM 

NSW 

NCAB 

NCUR 

ITER 

MTRC 

IBUOY 

ISPAR 

AXSM  (I) 

AYSM  (I) 

FSM  (I) 

FIDSM  (I) 

ASW  (J) 

FRSW  (J) 

FIDSW  (J) 

RHO 

AMC 

A  FAC 

TBH 

TBYMX 

SUBM 

TWX 

TIY 

TMIN 

CDASPX 

CDASBX 

VSP 

SPXK 

SPYK 

YYSP 

TIN  VI 

DTI 

TOTT 

DT2 

FLC  (K) 

DCI  (K) 

CDN  (K) 
CDT  (K) 

WC  (K) 

CM  (K) 
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APPENDIX  B  (Cont'd) 
CABUOY 

INPUT  PARAMETERS 
DATA  FORMAT  GUIDE 


11 

1-10 

FI  0.0 

Cl  (K) 

11-20 

FI  0.4 

C2  (K) 

21-30 

F10.4 

CINT  (K) 

31-40 

FI  0.2 

TREF  (K) 

41-50 

FI  0.2 

TENI (K) 

51-60 

FI  0.4 

PHID  (K) 

61-70 

FI  0.4 

XPI  (K) 

71-80 

FI  0.4 

YPI  (K) 

12 

1-10 

FI  0.4 

WBD  (K) 

11-20 

FI  0.4 

CDABX  (K) 

21-30 

FI  0.4 

CDABY  (K) 

31-40 

FI  0.4 

XMBV  (K) 

41-50 

FI  0.4 

YMBV  (K) 

13 

1-10,... 

FI  0.4 

YY  (L) 

14 

1-10,... 

FI  0.4 

CCK  (L) 

15 

1-10 

FI  0.4 

CDASY 

11-20 

FI  0.4 

WAST 

21-30 

FI  0.4 

RWY 

31-40 

FI  0.4 

RTX 

41-50 

FI  0.4 

RTY 

51-60 

FI  0.4 

YCG 

61-70 

FI  0.4 

BINT 

16 

1-10 

FI  0.4 

XSI 

11-20 

F10.4 

ZETI 

21-30 

F10.4 

SYDI 

31-40 

F10.4 

XPSI 

41-50 

FI  0.4 

ZTPI 

51-60 

FI  0.4 

SYPDI 

61-70 

FI  0.4 

DFTLIM 

71-80 

This  card  is  only  needed  if  2900<  FSM  (1 )  <  3000 

FI  0.4 

SYDLIM 

17 

1-10 

FI  0.4 

AKZ 

11-20 

FI  0.4 

AXP 

21-30 

FI  0.4 

APP 

31-40 

FI  0.4 

AFKX 

41-50 

FI  0.4 

AFKZ 

51-60 

FI  0.4 

AFKP 

B-2 
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APPENDIX  C 

CABUOY  OUTPUT  VARIABLES 


Ocean  Conditions 
Surface  Wave 

FREQ  —  Frequency  of  wave  component  (hz) 

AMPL  —  Amplitude  of  wave  component  (feet) 

PHASE  —  Phase  of  wave  component  (degrees) 

WL  —  Wave  length  of  component  (feet) 

WK  —  Wave  number  of  component  (1/feet) 

Current  Profile 

DEPTH  —  Depth  of  current  data  point  (feet) 

CURR  —  Current  velocity  at  depth  (DEPTH)  (knots) 

Surface  Motion 

FREQ  <  1000 

FREQ  —  Frequency  of  surface  motion  component  (hz) 

X-A  —  X  amplitude  of  motion  component  (ft) 

Y-A  —  Y  amplitude  of  motion  component  (ft) 

PHASE  —  Phase  shift  of  motion  component  (both  X&Y)  (degrees) 

2000  <  FREQ  <  3000 

FREQ 
X-A 
Y-A 
PHASE 

3000<  FREQ 

FREQ  —  Input  data  for  FSM  (Spheroid) 

X-A  —  Radius  of  float  at  waterline  (ft)  in  X-Z  plane. 

Y-A  —  Radius  of  float  at  waterline  (ft)  in  Y-Z  plane. 

PHASE  —  0 

General  Cable  Characteristics 

AM  COEFF  —  Cable  Added  Mass  Coefficient  (AMC) 

AREA  FAC  —  Cross  Sectional  Area  Correction  Factor  (AFAC) 

T  MIN  —  Minimum  tension  cable  can  support  (lbs) 


—  Input  numbers  for  FSM  -  Spar  buoy 

—  Cross  sectional  area  at  distance  Y-A  (ft  *)  from  origin 

—  Position  of  cross  sectional  area  data  point  relative  to  origin  (ft). 

—  0 


C-1 


NADC-801  74-30 


APPENDIX  C  (Cont'd) 


Cable  Properties 


Body  Properties 


NUM  —  Cabie  Segment  Number 

LEN  —  Length  of  segment  at  reference  tension  (ft) 

DIAM  —  Diameter  of  segment  (in.) 

CDN  —  Coefficient  of  drag  normal  to  cable 
CDT  —  Coefficient  of  drag  tangent  to  cable 
W  —  Wet  weight  of  cable  (Ibs/ft) 

M  —  Mass  of  cable  (slugs/ft) 

T  —  Reference  tension  of  cable  (lbs) 

Cl.  EXPC2  Coefficients  for  stress  strain  equation  T  =  clc 
CINT  —  Internal  damping  coefficient 


C2 


CDAY 

CDAX 

WT 

XV  M 

YVM 

X-Load 

Surface  Buoy  Characteristics 

WIND  LD 
SUBM 
CDASPX 
CDASBX 

AVSPX  1 

AVSPY  ) 

MAXTEN 

BUOY  CALC 

SPAR  NUM 

CDAY 

WT 

RWY 

RTX 

RTY 

YCG 

IN 

XI 

ZETA1 
SYI  ) 

XVI  * 

ZTVI  1 
SYVI  I 


—  Drag  area  in  X-Z  plane  (ft^) 

—  Drag  area  in  Y-Z  plane  (ft^j 

—  Wet  weight  of  body  (lbs) 

—  Virtual  mass  in  X  direction  (slugs) 

—  Virtual  mass  in  Y  direction  (slugs) 

on  BOT-WT  =  Applied  force  on  lower  weight  (lbs)  (TBH) 


—  Wind  loading  on  surface  float  (lbs)  (TWX) 

—  Depth  at  which  current  acts  on  surface  float  (ft) 

—  Drag  area  of  subsurface  package  in  x  direction  (ft  *) 

—  Drag  area  of  surface  float  excluding  subsurface  package  in  x 
direction  (ft^) 

Virtual  mass  of  surface  package  in  x  and  y  directions  (ft  ’) 

—  Maximum  tension  which  can  be  supported  by  the  cable  (lbs) 

—  variable  IBUOY 

—  Degree  of  polynomial  of  spar  flat  (ISPAR) 

—  Drag  area  in  X-Z  plane  (ft  2) 

—  Weight  in  air  of  float  (lbs) 

—  Vertical  distance  of  wind  loading  from  eg  (ft) 

—  Distance  in  x  direction  of  cable  attachment  point  from  cg  (ft) 

—  Distance  in  Y  direction  of  cable  attachment  point  from  cg  (ft) 

—  Submergence  of  cg  below  waterline  (ft) 

—  Moment  of  inertia  of  surface  float  (slugs  ft  squared) 

—  Initial  values  of  x.(.  ,<P  for  surface  float 
in  feet,  feet,  degrees 

—  Initial  velocities  of  x,  (  ,<t>  for  surface  float  in  feet/sec. 
feet/sec.,  degr  ;/sec 
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APPENDIX  C  (Corn'd) 


Constants  for  Spar  Buoy 

Initial  Ratio  TY  =  total  tension  at  top  of  cable  plus  weight  of  surface  float 
PO  =  p*  Surface  float  displaced  volume 
PI  =pjh  (y-y  Cg )  Sydy 

P2  =p/ £  (y-ycg)2  Sydy 
Q0  =p/^  (y-y  cg)  e*kV  Sydy 

Qi  =p/{]  (y-v  Cg)2e  ‘ky  svdv 
WVMZ  —  Virtual  mass  of  body  in  heave 
RGSO  —  p*g  ’cross  sectional  area  at  the  surface 
YCB-YCG  —  Separation  of  center  of  buoyancy  and  center  of  gravity 
EXACT  SYI  —  Steady  state  pitch  angle  of  float  in  degrees 
APPR  SYI  —  Value  of  float  pitch  used  in  dynamic  calculations.  May  vary 
from  exact  SYI  because  of  cable  segment  approximations 
made  in  CABUOY 

Constants  for  Spheroidal  Buoy 

DRAFT 
A 

VOL 
FKS 
FKH 
FKP 
FKPS 
RGSO 

CDAFC 

CDASX 
EXACT  SYI 
APPR  SYI 

Iteration 

ITN  —  Iteration  number 

VCXK  —  System  drift  velocity  (knots) 

TIXX  —  X  Component  of  tension  at  top  of  caole  (lbs) 

TIYY  —  Y  Component  of  tension  at  top  of  cable  (lbs) 

TENB  —  Tension  at  bottom  weight  (lbs) 

PHIB  —  Inclination  of  cable  at  bottom  relative  to  vertical  (degrees) 

TBXX  —  X  Component  of  tension  at  bottom  (lbs) 

TBYY  —  Y  Component  of  tension  at  bottom  (lbs) 


—  Draft  of  surface  float  (ft) 

—  Cross  sectional  area  at  the  waterline  (ft^) 

—  Displaced  volume  of  float  (ft  3) 

—  Added  mass  coefficient  in  surge 

—  Added  mass  coefficient  in  heave 

—  Added  mass  coefficient  in  pitch 

—  Coefficient  of  coupling  between  pitch  and  surge 

—  Water  density  *  acceleration  due  to  gravity  *  surface  float 
cross-sectional  area  at  the  waterline 

—  Drag  coefficient  for  submerged  portion  of  the  float  in  X 
direction 

—  Submerged  drag  area  of  surface  assembly 

—  Surface  fioat  tilt  (see  spar  buoy) 

—  Approximate  float  tilt  for  dynamic  calculations 
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APPENDIX  C  (Cont'd) 


XXBB  —  X  Coordinate  at  bottom  of  cable  (ft) 

YYBB  —  V  Coordinate  at  bottom  of  cable  (ft) 

SYSS  —  Pitch  of  surface  float  (degrees) 


Steady  State  Configuration 


TIX 

TIY 

DIRECTION 
NODE 
S  REF 
S  STR 
X 
Y 

TEN 

PHIS 


—  Tension  at  top  of  cable  in  X  direction  (lbs) 

—  Tension  at  top  of  cable  in  Y  direction  (lbs) 

—  Initial  conditions  known  at  top  or  bottom  (DIR) 

—  End  point  of  cable  segment 

—  Distance  along  reference  length  of  cable  (ft) 

—  Stretched  length  of  cable  (ft) 

—  X  position  of  NODE  (ft) 

—  Y  position  of  NODE  (ft) 

—  Tension  in  cable  at  NODE  (lbs) 

—  Tilt  of  cable  reference  to  vertical  (degrees) 


Cable  Initial  Conditions 


PHI 

XV 

YV 


—  Tilt  of  cable  at  NODE  (degrees) 

—  X  Velocity  of  NODE  (ft/sec) 

—  Y  Velocity  of  NODE  (ft/sec) 


Computer  Cable  System  Motions 


T 

DT 

X 

Y 

XP 

YP 

XPP 

YPP 

TEN 

FI 

FIP 

STRAIN 

STP 

SYPP 

DF 

PCT  T  DFTLIM  \- 
PCT  T  SYDLIM j 
WAVE 
BUOY 
0 
1 
2 


Current  time  of  dynamic  data  (sec) 
Integration  time  step  (sec) 

X  Position  of  body  (ft) 

Y  Position  of  body  (ft) 

X  Velocity  of  NODE  (ft/sec) 

Y  Velocity  of  NODE  (ft/sec) 

X  Acceleration  of  NODE  (ft/sec^) 

Y  Acceleration  of  NODE  (ft/sec^) 

Cable  Tension  (lbs) 

Angle  of  cable  at  NODE  (degrees) 

Rate  of  change  of  angle  (deg/sec) 

Cable  strain 

Rate  of  cable  strain  (1/sec) 

Angular  acceleration  of  float  pitch  (eg/sec^) 

■  Surface  float  draft  (ft) 

•  Percent  of  time  float  exceeds  draft  and 
Tilt  limits 

■  Data  for  surface  wave 

■  Data  for  origin  of  surface  float 

■  Data  for  top  of  cable  segment  number  1 

•  Data  for  bottom  of  cabie  segment  number  1 

■  Data  for  bottom  of  cable  segment  number  2 


Note:  Data  for  bottom  of  last  cable  segment  is  not  printed  out  as  point  is  assumed  to  be  fixed  in  space 
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